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'The Lord of the rings, the One Ring to rule them all' 
from J. R. R. Tolkien (D 

The standard theoretical basis for understanding superfluidity in Bose systems was formulated by Beliaev in 
1957, based on splitting the quantum field operator into a macroscopically occupied condensate component and 
a non-condensate component. This leads to a description of the condensate in terms of a 'single-particle state', 
the so-called macroscopic wavefunction. Since the discovery of Bose-condensed gases, an alternative theoretical 
picture has been developed which is based on a 'coherent band' of classically occupied states. This is often 
called the classical or c-field approach. The goal of this chapter is to review the differences between the Beliaev 
broken symmetry and c-field approach, and to argue that the c-field concept of a coherent condensate band of 
states has problems as a description of Bose superfluidity. However, the c-field idea of treating the lowest energy 
excitations classically can be used to advantage to simplify calculations within the Beliaev broken-symmetry 
formalism. 



I. INTRODUCTION 

One of the major goals of the Sandbjerg and the Durham 
workshops was to bring together theorists who were using dif- 
ferent methods for dealing with Bose gases at finite tempera- 
tures and in non-equilibrium states. The hope was that this 
would lead to a critical appraisal of the different formalisms, 
possibly including 'benchmark' problems where the results 
could be compared. These goals have been partially met. In 
particular, Proukakis and Jackson |2| have published a very 
useful article comparing many different formalisms for treat- 
ing Bose superfluids at finite temperatures. However a major 
problem in the recent literature is that papers based on the 
ZNG formalism |[3] 13 and the classical field approach |l5}d 
are like 'two ships passing in the night'. Apart from a few 
references, papers using one approach largely ignore results 
using the other approach. A major review article |6| on the c- 
field approach has 223 references, but the key paper setting up 
the ZNG formalism |3| is mentioned in one paragraph. By the 
same token, the recent monograph |4J by Griffin, Nikuni and 
Zaremba (henceforth refeiTed to as GNZ) only makes brief 
contact with the final results of c-field theories in Chapters 8 
and 9 when deriving the simple dissipative Gross-Pitaevskii 
Equation (DGPE). 

In an attempt to remedy this situation and gain greater un- 
derstanding, we initiate in this chapter a critical discussion of 
the conceptual basis of the classical field approach |l2l|5]|6l[8l, 
comparing it with the Beliaev broken-symmetry theory of 
superfluids (see Ref. |9| and also Chapters 4-5 of GNZ). 
The latter remains the basis of almost all current discus- 
sions of superfluidity in interacting Bose systems in the lit- 
erature using quantum field theoretic techniques (see, for ex- 
ample, Ref. |10|). The ZNG coupled equations [Eqs. (3.21) 
and (3.42) of GNZ] are based on the Beliaev scenario, albeit 
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with the addition of various approximations ('Popov' approx- 
imation for the neglect of the pair anomalous average fTT), 
Hartree-Fock spectrum for the excitations, etc.). However, as 
discussed in Chapters 6 and 7 of GNZ, the simplified ZNG 
formalism can be extended in a systematic manner within the 
Beliaev field-theoretic approach. 

We hope our comparison of the basis of the c-field approach 
with the Beliaev formulation will make a useful contribution 
to the current literature. For references and more detailed dis- 
cussion of the topics covered in this chapter, we refer to recent 
reviews [2^ 4^-61 . An alternative viewpoint on this important 
comparison is given in Ref. lfT2l . 

II. METHODOLOGY 
A. Beliaev Theory of Superfluidity 

The key feature of any theory of a Bose-condensed quan- 
tum fluid is that it should explain the appearance of a new 
degree of freedom which describes the phenomenon of super- 
fluidity. These features were first captured in the Landau two- 
fluid theory (developed in 1941 for superfluid "^He) and were 
later given a full microscopic basis by Beliaev in 1957 |9|. Be- 
liaev, extending the pioneering work of London fT3l in 1938 
and Bogoliubov [14J in 1947 on dilute Bose gases, formulated 
a systematic way of separating out the superfluid degree of 
freedom associated with the appearance of a Bose condensate 
(BEC) in a quantum field-theoretic description of an interact- 
ing system of Bose particles. In Beliaev's formulation, which 
is the canonical approach to superfluidity in Bose systems, the 
superfluid degree of freedom has its microscopic origin in the 
broken-symmetry average value of the quantum field opera- 
tor 0(r) = (*I'(r))BS- This can be viewed as a single-particle 
state which is macroscopically occupied and plays the role 
of the order parameter of the new superfluid phase appearing 
below some transition temperature T^. The concept can also 
be extended to non-equilibrium situations where the ensem- 
ble average is taken with respect to a non-equilibrium density 
matrix. 
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The Beliaev formulation f9] is based on a decomposition of 
the quantum field operator *J'(r) into a condensate component 
0(r) as defined above plus a non-condensate component'. 



Y(r) = 0(r) + ,A'(r). 



(1) 



Invoking this idea, one can then develop a many-body per- 
turbation scheme for the non-condensate fluctuation (^'(r) in 
the presence of the Bose-condensate described by 0(r). This 
takes into account fully the complex correlations between 0(r) 
and (i'Xr). Because it describes a single-particle state which is 
macroscopically occupied, in most treatments (p(r) is treated 
as a classical field. However this latter approximation is not 
essential to the Beliaev broken-symmetry formalism llT5irT6l . 

A crucial aspect of the Beliaev formulation of superfluidity 
is that the macroscopic wavefunction 0(r) can be complex 



0(r) = ^/^)e''>^'\ 



(2) 



where ndr) - |0(r)p defines the condensate density. The 
phase of this order parameter naturally leads to an associated 
condensate velocity field 



Vc(r) = -V0(r). 

m 



(3) 



Clearly Vc(r) is irrotational, a key element of the Landau the- 
ory of superfluidity. Choosing a finite value of the order pa- 
rameter 0(r) corresponds to a specific choice of the phase, 
which breaks the gauge symmetry. This is most conveniently 
done by using a number non-conserving formalism 14] |9] \T5\ - 

Using linear response theory to describe the effect of a con- 
densate moving with a velocity Vc, one can calculate how 
it drags the non-condensate atoms along. The total super- 
fluid mass current can be rigorously proved [see Eq. (3.3) of 
Ref [ 15)] to be given by js = PsVc, where the superfluid den- 
sity Ps(r) is fundamentally distinct from ndr) and is given as 
the difference between the longitudinal and transverse veloc- 
ity response functions (seepages 128-131 of Ref. 1171 ). Thus 
the Beliaev broken-symmetry formalism based on Eqs. ([T]i 
and (j2]i naturally leads to two fundamental features of a Bose 
superfluid, namely (a) the superfluid velocity Vs is equal to the 
condensate velocity Vc, and (b) a rigorous microscopic defi- 
nition for the superfluid density, which reduces to Landau's 
quasiparticle expression in two-fluid hydrodynamics. These 
two features do not have to be inserted into the theory as ad- 
ditional assumptions. 

The Beliaev picture, based on Eq. ([T]), is the conceptual and 
computational basis of most microscopic studies of superflu- 
idity in liquid "^He. The work of ZNG [4 1 as applied to trapped 
Bose gases is based on this approach, albeit within a simpli- 
fied treatment of the non-condensate fluctuations. One of the 
main strengths of the Beliaev approach is that it provides a mi- 



croscopic basis for the derivation of Landau's phenomenolog- 
ical two-fluid hydrodynamics which is valid when interactions 
are strong enough to establish local thermodynamic equilib- 
rium. Likewise, it can be used to derive Landau-Khalatnikov 
hydrodynamics (see Chapter 17 of GNZ) which accounts for 
the damping that is associated with dissipative transport co- 
efficients. In this derivation, a natural definition of what we 
'mean' by the superfluid density Ps(r) emerges, based on the 
underlying Bose condensate density «c(r). While Ps(r) can be 
quite different from «c(r), as noted above, the superfluid ve- 
locity field Vs(r) is equal to the condensate velocity field Vc(r) 
defined in Eq. ([3]l. Being able to derive Landau's two-fluid 
equations under appropriate conditions (short collision times) 
is an essential requirement of any complete microscopic the- 
ory of superfluidity. 

It is clear that in the Beliaev formulation, there is a unique 
single-particle state described by 0(r) which has a 'privileged 
status' among other single-particle states and describes the 
new superfluid phase. To quote J. R. R. Tolkien when he refers 
to the 19 rings produced and distributed by Sauron of Mor- 
dor HI, one might say that 0(r) is 'The Lord of the rings, the 
One Ring to rule them all'. In recent years, however, a seem- 
ingly difi'erent picture of describing interacting Bose gases has 
emerged 111 |5] ID [8^ in which 0(r) has no special status. This 
approach is now often called the 'classical-field approach', or 
more succinctly, the c-field approach. In Section II B we re- 



view the conceptual basis of this approach. In Sections III A 
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and |IIIB| we then argue that despite its success as a compu- 
tational method, it has problems when it comes to addressing 
superfluidity in interacting Bose systems. 



B. C-Field Approach 

Recently there has been considerable work on Bose gases 
at finite temperatures using what is called the c-field approach 
(see Refs. ifTSl and ||T9| for a general overview). The key idea 
behind this approach is based on the observation that the low- 
est energy eigenstates of a trapped Bose gas have occupation 
numbers much larger than unity and hence can be treated clas- 
sically fSG'-^Tl. One variant of this approach f6\ divides the 
states of an interacting Bose gas into two classes: a 'coherent' 
band of states which is treated by means of a classical field, 
and an 'incoherent' band of higher energy states which must 
be treated in some other way. The classical field can be intro- 
duced systematically by means of a projection technique (as 
discussed further below). The c-field approach is an attractive 
numerical scheme since the states in the coherent band satisfy 
an equation analogous to the usual Gross-Pitaevskii Equa- 
tion (GPE) valid at T = 0. As noted above, a very detailed 
review article ||6l has been published discussing the c-field 
approach. It sets up the conceptual basis for this approach, 
with extensive references to earlier work, and shows how var- 
ious kinds of projected and stochastic Gross-Pitaevskii-type 
equations can be derived. The projected Gross-Pitaevskii 
equation (PGPE) omits entirely the coupling to the incoher- 
ent region |24|, whereas the so-called stochastic theories in- 
clude it via dissipation and stochastic noise terms. These 
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latter theories include the closely-related stochastic projected 
GPE (SPGPE) of Gardiner et al. lH and the stochastic GPE 
(SGPE) approach of Stoof ||26l . In the following, we focus 
on the conceptual issues underlying the PGPE approach (see 
Refs. L5J and [6J for details regarding the numerical imple- 
mentation of the c-field scenario), although many of our com- 
ments are also relevant to the SPGPE and SGPE approaches. 
We refer separately to the SGPE approach of Stoof in Sec- 
tion |IlID] 



In Section |II A| we pointed out that the Beliaev formula- 
tion for dealing with a Bose system of particles with a con- 
densate is based on the separation of the quantum field op- 
erator 4'(r) into two distinct parts, as given in Eq. The 
end result is a rigorous scheme for treating the dynamics of 
the non-condensate field operator (^'(r) which is dynamically 
coupled to the superfluid degree of freedom <pir). This opera- 
tor of course also describes the non-condensate in the normal 
phase. The c-field approach adopts a somewhat similar point 
of view, except that the quantum field operator *I'(r) is now 
decomposed in a manner different from Beliaev, namely 



Y(r) = ^c(r) + ^PiCr). 



(4) 



Here the part of which is referred to as the coherent part 
4*0 corresponds to a band of single-particle states which can 
be treated classically since their occupation numbers are as- 
sumed to be much larger than unity. The incoherent part 
4*1 describes the higher energy states which must be treated 
quantum mechanically. The best choice of 'energy cutoff' 
Ecut which separates these two bands is a delicate question 
discussed at length in the implementation of the c-field ap- 
proach [2, 5 , 6|, but one might expect E^ut ~ ksT- 

The simplest (and earliest |'20'-'23l) version of the c-field ap- 
proach is to restrict oneself to the coherent 'condensate' band. 
The exact Heisenberg equation of motion for *I'(r) is given by 



ih 



d'¥{r, t) 
dt 



h^v^ 

2m 



+ Vext(r) ^^(r, t)+g^Hr, tyVir, ffVir, t), 

(5) 

where we assume a simple i-wave delta function pseudopo- 
tential of strength g. Writing *I'(r) as an expansion in an ap- 
propriate basis of single-particle eigenstates </J,(r), 



4'(r) = 2^,(r)fl, 



(6) 



the coherent part in Eq. Q can be expressed as the projection 
onto the states comprising the coherent band 

*c(r) = !Pc{*(r)} 

= J]<Pi(r) [ «'rV;(r'mr'). (7) 

To the extent that this coherent part of the full quantum field 
operator describes states which can be treated classically [i.e., 
4'c(r) becomes a c-number function], it satisfies the equation 



of motion [cf Eq. Q] 



ih- 



dt 



+ Vex,(r) 



2m 

{g|*I'c(r,f)l'^c(r,f) 



(8) 



This projected Gross-Pitaevskii equation (PGPE) is quite dif- 
ferent from the usual GPE since *I'c(r, describes a band of 
states. However, the techniques for solving Eq. ([8]) are very 
similar to those used to solve the GPE for 0(r, f) at T = 0. 

An essential difference between the c-field approach and 
the Beliaev formulation is already apparent. *I'c(r) describes 
a 'band' of low energy single particle eigenstates, rather than 
the lowest one ^{v). This band is sometimes referred to as 
a 'coherent' or 'condensate' band of states and plays a cru- 
cial conceptual role in c-field treatments. Equally important, 
it forms the basis for numerical approximations [2 , 5 , 6 1. It is 
correctly argued that this condensate band (even when the in- 
coherent part is ignored) contains an essential aspect of the 
physics of Bose superfluids at finite temperatures which is 
omitted when only the average 0(r) is considered, namely the 
fluctuations of the condensate. 



III. VALIDITY 

A. Conceptual Problems with the C-Field Approach 

In our view, a major shortcoming of the c-field theories is 
that no fundamental distinction is made between the different 
low energy states contributing to the coherent band. Indeed, 
the classically occupied states comprising the coherent band 
are sometimes viewed as defining a 'kind of order parame- 
ter' which generalises the single state ^(v) used in the Beli- 
aev broken-symmetry theory (see, for example, Refs. [2,, ^). 
In such treatments, it is not made clear how this band of 
states can be considered an 'order parameter' in the same 
way that ^{r) in the symmetry-breaking scenario is associ- 
ated with a new superfluid phase with off-diagonal long-range 
order (ODLRO). It simply picks out a class of low energy 
states which can be treated classically. This set of states is 
certainly not the 'condensate' in the conventional sense (as 
also highlighted in those works). In order to make contact 
with the condensate in the Bose broken-symmetry picture, the 
c-field methods invoke a Penrose-Onsager analysis Il27l l28ll 
to diagonalise the single-particle density matrix p^''(r, r', t) - 
("i*' (r, t)'^{r' , t)) where the average is to be taken with respect 
to a non-equilibrium density matrix. In the c-field (PGPE) ap- 
proach, *i'(r, i) is replaced by *I'c(r, and the average is eval- 
uated as a time average. The lowest eigenstate is then taken to 
define the Bose condensate. This approach is clearly reason- 
able for obtaining the condensate in the equilibrium state but 
some different method must be used to make contact with the 
time-dependent order parameter 0(r, t) that would appear in a 
non-equilibrium situation. The latter quantity appears directly 
when Bose broken-symmetry is invoked in non-equilibrium 
theories such as the ZNG theory. However, it should be noted 
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that in the context of the SGPE, 0(r, t) can also be extracted by 
performing an additional ensemble average over a distribution 
of initial states consistent with the non-equilibrium evolution 
of interest (see Section [riID[ ). 

As emphasised in our review of the Beliaev formulation in 
Section II A the gradient of the phase of <p{r) introduces a 



unique superfluid velocity field Vs(r). The existence of such a 
velocity (sometimes viewed as a 'phase locking') defines what 
we mean by superfluidity and is the reason Beliaev's formu- 
lation is useful in the microscopic derivation of the two-fluid 
behaviour of Bose-condensed systems. In contrast, each of 
the states arising from the Penrose-Onsager diagonalisation 
of the 'coherent band' would in general be expected to have a 
different velocity field, related to the phase of each eigenstate. 
It of course would be natural to identify the superfluid velocity 
with the macroscopically occupied state, that is, the one corre- 
sponding to the lowest Penrose-Onsager eigenstate. However, 
in recent expositions 121 15] |6| of c-field theories there is little 
discussion of the superfluid velocity field, a concept central 
to any satisfactory theory of a Bose superfluid. Whenever the 
superfluid velocity is addressed there necessarily is a 'Lord of 
the rings'. 



B. Excitations in the C-Field and Beliaev Approaclies 

It should be noted that the power of Beliaev's 'extraction' 
of the condensate part of 0(r) as defined in Eq. ([T]i lies in the 
fact that one is able to set up a systematic diagrammatic the- 
ory for the non-condensate dynamics completely analogous 
to that used for a Bose gas in the normal phase @[l6l|T71. It 
is far from obvious that a similar perturbation scheme for the 
incoherent band *I'i(r) can be developed based on the c-field 
decomposition given in Eq. Q. 

Reading the literature on Bose superfluids, one can eas- 
ily get the wrong impression that the excitation spectrum of 
the condensate and non-condensate components [obtained via 
Eq. ([T]l] are different. In fact, one of the great triumphs 
of the field-theoretic analysis of the structure of correlation 
functions is that it shows explicitly how a condensate cou- 
ples and hybridises the excitations of the condensate micro- 
scopic wavefunction 0(r, t) with those of the non-condensate 
described by the single-particle Beliaev Green's functions for 
(/('(r). In the field-theoretic formalism (see pages 69-74 of 
GNZ), one has for example the decomposition 

Gi(i,i') = -/<r^(i)^"'(i')> 

= -/<«A'(i)-A'^d'))+ V^0(i)V^0*(i') 

= Gi(1,1') + Gi/2(1)G;/2(1'). (9) 

Bose broken-symmetry leads to equations of motion for Gi 
and Gi/2 which are coupled and this in turn leads to both of 
these functions sharing the same hybridised excitation spec- 
trum. In an analogous manner, the single-particle Green's 
function Gi(l, 1') and the density response function ;t'„„(l, 1') 
share the same excitation spectrum (see pages 96-98 of GNZ), 
a key experimental signature of Bose superfluids. This com- 



mon excitation spectrum (induced by the condensate) is im- 
plicitly contained in the pioneering Green's function analysis 
of Hohenberg and Martin (see Section VI of Ref. ||29l ) and 
was exhibited more explicitly at finite temperatures by Che- 
ung and Griffin |,30J . A more general way of showing this 
shared excitation spectrum is to use the diagrammatic dielec- 
tric formalism, reviewed in Chapter 5 of Ref. fTTl. For a 
model calculation showing how the hybridising eff'ect of the 
condensate couples the excitations of the condensate and non- 
condensate, see Section 5.4 of GNZ. 

This shared excitation spectrum is usually hidden in most 
theoretical papers since G1/2, Gi, andxm are computed at dif- 
ferent levels of approximation. As a result, the excitations 
associated with these fluctuations appear to be different. In 
a consistent calculation, as discussed above, all these corre- 
lation functions exhibit the same hybridised excitation spec- 
trum — the characteristic signature of a Bose superfluid. To 
be explicit, the GPE-like equation of motion for the conden- 
sate Gi/2 (i.e., (p) will have the same excitation spectrum as 
the non-condensate correlation function Gi, but only if we 
take care to use consistent approximations in both equations 
of motion. 

The preceding discussion indicates that the c-field division 
of the excitations of a Bose superfluid into two regions, the co- 
herent and incoherent bands, is both artificial and misleading. 
In the Beliaev formulation, there is a single excitation branch 
describing both the condensate and non-condensate fluctua- 
tions. At low momentum, these excitations in a Bose super- 
fluid always have the characteristic spectrum of Goldstone- 
Nambu phonons. At high momentum, the excitations are 
particle-like in a dilute Bose gas (in contrast to rotons in su- 
perfluid "*He) and can be described using a Hartree-Fock (HE) 
spectrum. Which momentum region has the greatest weight in 
thermodynamic quantities depends on whether one is consid- 
ering the very low (phonon-like excitations) or high (particle- 
like excitations) temperature region of the superfluid phase. 

In this regard, the ZNG model calculations carried out to 
date treat the fluctuations of lA'(r) in terms of HE-like exci- 
tations moving in a time-dependent self-consistent field pro- 
duced by the condensate and non-condensate components. 
Although the condensate degree of freedom and the thermal 
cloud of atoms are treated using quite different approxima- 
tions, the ZNG theory nevertheless gives an excellent account 
of the frequencies and damping of collective modes over the 
range of temperatures studied experimentally (see Chapters 1 1 
and 12 of GNZ for a review of applications). The essential 
reason for this is that the dynamics of the condensate and 
non-condensate are coupled in a physically realistic way. As 
discussed in Chapter 7 of GNZ, one can generalise the ZNG 
equations of motion by treating the thermal cloud atoms in 
terms of Goldstone-Nambu phonons, thereby extending the 
validity of the theory to lower temperatures and restoring the 
constraints on the excitation spectrum imposed by Bose bro- 
ken symmetry. However, from a practical point of view, such 
an extension would have little effect on the calculated damp- 
ing rates at higher temperatures. 

In the review article |2| comparing different methods for 
dealing with Bose gases at finite temperatures, it is suggested 
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that the Behaev formulation (on which ZNG is based) rests on 
an 'artificial' separation into condensate and non-condensate 
contributions (the latter being the thermal cloud atoms in 
ZNG). We maintain that this separation is not artificial but 
gives a natural way of capturing the two-fluid nature of su- 
perfluids resulting from an underlying Bose condensate. In 



contrast, as discussed here and in Section III A the separation 
of the states into coherent and incoherent bands in the c-field 
approach is somewhat ad hoc. It leads neither to an under- 
standing of the special role of the lowest energy condensate 
mode and the superfluid velocity field associated with it, nor 
of the special nature of excitations in a Bose-condensed fluid. 

A further limitation of the c-field theories when including 
the coupling to the incoherent band of high energy states (i.e., 
in the context of the SPGPE) is that the latter is approxi- 
mated as an effective heat bath in all (numerical) treatments 
to date, while the microscopic dynamics of the coherent band 
is worked out in detail. Without a proper treatment of the dy- 
namic coupling between the coherent and incoherent bands 
it will be impossible to describe the damping of collective 
modes in an accurate way. As stated earlier, the natural way to 
include this coupling is by means of a many-body perturbation 
theory approach based on Bose broken symmetry. 



C. Computational Advantages of the C-Field Approach 

Despite our critical remarks regarding the conceptual ba- 
sis of the c-field approach, we must concede that being able 
to treat the low-lying states classically is an important com- 
putational advance. As a result, one can use a PGPE [such 
as Eq. ([8])] to address some non-trivial non-equilibrium prob- 
lems without having to introduce a lot of formal 'machinery' 
typical of field theoretic calculations based on the Beliaev 
formalism (diagrammatic perturbation theory, single-particle 
Green's functions, etc.). By its nature, the c-field method is 
able to provide information regarding fluctuations not readily 
available using other approaches. In addition, to the extent 
that the incoherent band can be treated as a thermal bath as 
done in the SPGPE, it can provide a qualitative understanding 
of various equilibration processes such as vortex nucleation 
and relaxation in rotating condensed gases. It would never- 
theless be useful to have more detailed comparisons with other 
approaches for situations where the dynamics of the incoher- 
ent band cannot be neglected. 

Another interesting question is to what extent the c-field 
approach can provide an adequate description of the critical 
region near T^. It is well known from the theory of second 
order phase transitions (arising from a broken-symmetry or- 
der parameter) that fluctuations not included in a mean field 
approximation are crucial. This topic is surprisingly difficult 
and requires a careful treatment of the infrared divergences 
that arise when using quantum field theoretic methods (see, 
for example Refs. IfTOl [3T]). In the c-field literature, it has 
been argued that the PGPE for the coherent band naturally 
includes the low energy critical fluctuations near (see, for 
example, Refs. I6l [32l ). At some level, this conjecture is rea- 
sonable. However, one would like more analytical studies to 



determine if such PGPE calculations (which have the virtue of 
simplicity) indeed capture the results found in the more stan- 
dard many-body literature IfTOl [311 . 

In summary, we believe the c-field approach as currently 
implemented is best viewed as a numerical strategy to sim- 
plify calculations rather than as an 'alternative' theory to the 
Beliaev approach. 



D. Stoof Formalism: Some Brief Remarks 

Apart from ZNG and the c-field approach, the only 
other general scheme which attempts to deal with the non- 
equilibrium properties of a Bose-condensed gas is the path- 
integral formalism developed by Stoof |l8]|33l. The final equa- 
tions he derives are closely related to those of the ZNG ap- 
proach. However, by formulating a theory for the probability 
distribution of a bosonic field, fluctuations of the condensate 
are included. A simplified version of his theory results in a 
stochastic GPE (SGPE) with a noise term The pres- 

ence of the noise term distinguishes this SGPE from the PGPE 
of the c-field method discussed earlier; in fact, the more recent 
formulations of the c-field method (SPGPE) are essentially 
equivalent to Stoof's theory (see Ref. |2|). 

In applying this general formalism, however, it is important 
to note that Stoof and coworkers make two further simplifi- 
cations. As in the ZNG theory, the thermal cloud excitations 
are treated within a simple self-consistent Hartree-Fock ap- 
proximation. However, more importantly (for the purposes 
of numerical implementation) the non-condensate component 
(thermal cloud) is assumed to be in static thermal equilib- 
rium, described by the usual equilibrium Bose-Einstein dis- 
tribution. As noted earlier, this assumption is also standard in 
current applications of the c-field approach which include the 
coupling to the incoherent region (SPGPE). However, explicit 
calculations using the ZNG coupled equations (see Chapter 12 
of GNZ) have demonstrated that the inclusion of the thermal 
cloud dynamics is essential in the determination of the fre- 
quency and damping of condensate collective modes at higher 
temperatures. The static thermal cloud approximation (see 
Chapter 8 of GNZ) may be useful as a starting point in under- 
standing the eff'ect of the thermal cloud, but it does not lead to 
quantitative predictions. 

One might argue that the neglect of the dynamic HE mean 
field arising from the non-condensate atoms is a reasonable 
approximation since it is small in comparison to the conden- 
sate mean field. However, we emphasise that the neglect of the 
thermal cloud HE mean-field 2gn'{r, t) in the Stoof SGPE I8l 
is not a conceptual difference (as suggested in Ref. Q) with 
the generalised GPE of ZNG which includes it [see Eq. (3.21) 
of GNZ], but is simply an ad hoc simplification useful for 
carrying out explicit calculations within the Stoof formalism. 
By doing so, important physics is lost. For one example, as 
discussed in Chapter 13 of GNZ, this thermal cloud HF mean- 
field is the agent responsible for the Landau damping of con- 
densate collective modes. 
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